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ABSTRACT 

The stability of the random field Ising model (RFIM) against spin glass (SG) 
fluctuations, as investigated by Mezard and Young, is naturally expressed via 
Legendre transforms, stability being then associated with the non-negativeness 
of eigenvalues of the inverse of a generalized SG susceptibility matrix. It is found 
that the signal for the occurrence of the SG transition will manifest itself in free- 
energy fluctuations only, and not in the free energy itself. 

Eigenvalues of the inverse SG susceptibility matrix is then approached by 
the Rayleigh Ritz method which provides an upper bound. Coming from the 
paramagnetic phase on the Curie line, one is able to use a virial-like relationship 
generated by scaling the single unit length (D < 6; in higher dimension a new 
length sets in, the inverse momentum cut off). Instability towards a SG phase being 
probed on pairs of distinct replicas, it follows that, despite the repulsive coupling 
of the RFIM the effective pair coupling is attractive (at least for small values of the 
parameter gA, g the coupling and A the effective random field fluctuation). As a 
result, "bound states" associated with replica pairs (negative eigenvalues) provide 
the instability signature. Away from the Curie line, the attraction is damped out 
till the SG transition line is reached and paramagnetism restored. In D < 6, the 
SG transition always precedes the ferromagnetic one, thus the domain in dimension 
where standard dimensional reduction would apply (on the Curie line) shrinks to 
zero. 
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ABSTRACT 

The stability of the random field Ising model (RFIM) against spin glass (SG) 
fluctuations, as investigated by Mezard and Young, is naturally expressed via 
Legendre transforms, stability being then associated with the non-negativeness of 
eigenvalues of the inverse of a generalized SG susceptibility matrix. It is interesting 
to note that the signal for the occurrence of the SG transition will manifest itself in 
free-energy fluctuations only, i.e., in the n-replica formulation, in terms of order n 2 
and higher and not in the free energy itself (when the SG phase sets in, all terms of 
order n 2 or higher, become linear and start contributing to the free energy, a behavior 
contrasting with the standard SG). 

Eigenvalues of the inverse SG susceptibility matrix can be approached by 
the Rayleigh Ritz method which provides an upper bound. Coming from the 
paramagnetic phase on the Curie line, one is then able to use a virial-like relationship 
generated by scaling the single unit length (D < 6; in higher dimension a new length 
sets in, the inverse momentum cut off). Instability towards a SG phase being probed 
on pairs of distinct replicas, it follows that, despite the repulsive coupling of the RFIM 
the effective pair coupling is attractive (at least for small values of the parameter gA, g 
the coupling and A the effective random field fluctuation). As a result, "bound states" 
associated with replica pairs (negative eigenvalues) provide the instability signature. 
Away from the Curie line, the attraction is damped out till the SG transition line is 
reached and paramagnetism restored. In D < 6, the SG transition always precedes 
the ferromagnetic one, thus the domain in dimension where standard dimensional 
reduction would apply, shrinks to zero. 
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RESUME 

On examine la stabilite du modele d'Ising en champ aleatoire (MICA) par rapport 
aux fluctuations verre de spin (VS) comme chez Mezard et Young. Cette stabilite 
s'exprime naturellement via des transformations de Legendre ou elle est associee a 
la non-negativite des valeurs propres de la matrice inverse d'une susceptibilite VS 
generalisee. II est interessant de noter que le signal de la transition VS n'apparait que 
dans les fluctuations de l'energie libre, c'est-a-dire, dans la formulation en n-repliques, 
que dans les termes d'ordre n 2 ou plus eleve et non dans l'energie libre elle-meme (a 
la transition VS, les termes d'ordre n 2 ou plus eleve, se transforment tous en termes 
lineaire en n et contribuent alors a l'energie libre, un comportement qui contraste 
d'avec celui du VS standard). 

Les valeurs propres de la matrice inverse de la susceptibilite VS peuvent s'estimer 
par la methode de Rayleigh-Ritz qui fournit une borne superieure. Sur la ligne de 
Curie, en s'approchant a partir de la phase paramagnetique, on peut utiliser alors une 
relation du type "viriel" obtenue en dilatant 1' unique echelle des longueurs (D < 6; 
en dimension superieure une deuxieme echelle s'introduit, l'inverse du moment de 
coupure). L'instabilite VS se manifestant sur les correlations de paires de repliques 
distinctes il s'ensuit que, malgre le couplage repulsif du MICA, le couplage effectif 
d'une paire est attractif (du moins pour les faibles valeurs du parametre gA, ou g 
est le couplage et A la fluctuation effective du champ aleatoire). II en resulte des 
"etats lies" associes aux paires de repliques (valeurs propres negatives) qui sont la 
signature de l'instabilite. Quand on s'eloigne de la ligne de Curie, l'attraction de 
paires est amortie, jusqu'a la ligne de transition VS ou le paramagnetisme apparait. 
En dimension D < 6, la transition VS precede toujours la ferromagnetique, si bien 
que l'intervalle de dimension ou la reduction dimensionnelle pourrait s'appliquer, se 
reduit a zero. 
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After nearly twenty years of intense activity, there is yet no consensus on the 
critical behavior of random field systems (for recent reviews see [1,2]). 

A blatant contradiction arose when a calculation to all orders in perturbations! 3 ' 4 ] ,| 
later supported by a non perturbative approach! 5 ], established dimensional reduction 
(between the RFIM in dimension D and the pure Ising system in D — 9) both for hy- 
perscaling relationships between critical exponents and for the exponents themselves 
as a function of D. With 9 = 2 this was predicting a lower critical dimension Dg = 3 
for the existence of a ferromagnetic phase, in contradiction with an early Imry-Ma! 6 ] 
argument predicting Di = 2 (later supported by rigorous work of Imbrie^ , proving 
the existence of ferromagnetism in D = 3). 

Despite the fact that Dg = 2 is now widely accepted, there remains the question 
of down to which dimension are the resummed perturbation results valid, and what 
happens below that dimension. 

Meanwhile several groups have proposed the existence of a glassy phase sector 
in the A,T plane (A is the width of the random field gaussian distribution and T 
the temperature) out of numerical studies! 8-12 ] or from analytical work, extending to 
random field systems' 13-15 ] the techniques of replica symmetry breaking! 16 ] (RSB). 
In particular Mezard and Young! 14 ] have used as a control parameter the number of 
components m (as in Bray! 17 ] self consistent screening approximation) and written 
out explicit self consistent equations for the exponents. 

Here, we follow the most straightforward approach to analyze properties of the 
RFIM, i.e. we mimick what is done in the pure system to describe the paramagnetic 
and condensed (ferromagnetic) phases. 

We know since the work of Yvon! 18 ] that the appropriate way to have access to 
the condensed phase is to replace the expansion in the local field Hi, by one in the 
local magnetization Mj, through a Legendre transform. The Jacobian of the transform 
det (dMi/dHj) vanishes at the transition (with the lowest eigenvalue of the matrix 
dMi/dHj), displaying the non-equivalence of the Hi and the Mj expansions. 

Likewise here we consider the RFIM described by an effective hamiltonian with 
an external field A and perform the appropriate Legendre transform to the conjugate 
observable. Again the lowest eigenvalue of the jacobian matrix yields the locus of the 
singularities of the associated susceptibility, here the SG susceptibility, i.e. the line 
of the SG transition. 

In Section 1-2, we recall the perturbation expansion and effective hamiltonian 
for the RFIM. In Section 3, the Legendre transform is effected yielding stationarity 
conditions and eigenvalue equations for the SG transition. In Section 4 we study the 
transition and show that it manifests itself in the free-energy fluctuation (as contrasted 
with the standard SG). Section 5 is devoted to a study of the phase diagram using the 
Rayleigh-Ritz variational method! 19 ] and we conclude in Section 6, where our results 
are summarized. 
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1. THE PURE ISING SYSTEM 

Consider the pure Ising hamiltonian, in its soft spin version, 



« = \ E fa + p 2 ) <p(pM-p) + 1 E ^ - E tf^i (i-i) 



and 



{#,,} = ln J Dip exp - H{v?} = In Z (1.2) 
One may describe the system by expanding in H; or via a Legendre transform 

W {H 3 } = -r {Mj} + £ HjMj (1.3) 



where the magnetization Mj is 



9^ 

= = M ( L4 ) 



and the Legendre transform T satisfies 

by expanding in M. The bracket in (1.4) stands for "thermal" average 

Mj = J Dip ip 3 Q- n{ ^ jZ (1.6) 

The Jacobian of the transformation is the determinant of the inverse susceptibility 
matrix 

f^" 1 ) = (1.7) 

and when a zero eigenvalue occurs it signals the inequivalence of the two expansions 
and the occurrence of a transition. In Fourier transform x _1 ((/), and e.g., in zero 
momentum for a standard system, x _1 vanishes in zero field, at T c the Curie point 

X - 1 (q = 0;T c ) = 0, (1.8) 

below which M^O even for H = 0. 

This well known description of the paramagnetic to ferromagnetic transition we 
would like now to extend to the random field system. 
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2. THE RANDOM FIELD ISING SYSTEM: Perturbation, a reminder! 20 ] 

Let us consider now Hi to be a quenched random field with a pure gaussian 
probability distribution i.e. 

Ht=0 (2.1) 



H i H j = SijA (2.2) 
where the bar stands for probability average. 

2.1 Direct averaging 

One may compute the H expansion of observables and then perform the Wick 
average of the if's on each term of the expansion. In the paramagnetic phase, to keep 
things simple, one obtains 



W{H} = ^ connected graphs with all pairs of H's 

coalesced as in (2.2). (2.3) 



G(i;j) = (<Pi<Pj) — ((pi) (<Pj) = connected graphs, rooted at i and j, 

with all pairs of H's coalesced. (2.4) 



C(i;j) = (cpi) (ifj) = graphs made of two disconnected pieces, respectively rooted 

at i and j with all pairs of H's colaesced. (2.5) 

The coalescence of all pairs of H } s transforms the two disconnected pieces into 
one single field-connected graph. 

2.2 Averaging via the effective hamiltonian 

The above results can be recovered using the replica trick, that is, computing 



2 

rr 



Z n = (exp W{H}) n = exp { n W + — 



where one then recovers the averaged free energy 



-F = W = [ Z n - 1 



, \ 2 



w 2 -yw 



+ ... (2.6) 



(2.7) 
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but also its successive fluctuation cumulants. 
The effective hamiltonian is now, 

"» = E ( 5 E Co + f 2 ) *> w(-p) + 1 E (^°) 4 } - 5 E E A *> 

a ^ p j J a,/3 p 

(2.8) 

The propagator becomes a matrix G with components 

G a( 3 = {<Pa<Pp) n - (<f a ) n (<^3> n (2-9) 

where the replica-thermal average is shown as () n . In the paramagnetic phase (no 
magnetization) ((f a ) n = 0. The bare propagator is the inverse of the matrix 

(p 2 + t ) 5 a p - A, i.e. with [G°] _1 = p 2 + t 

g %- g "^ + T^g«- (2 ' 10) 

The first term is the connected (bare) propagator, the last is the field-connected (bare) 
propagator (suppressing the H- coalescence into A's it falls into several disconnected 
pieces). In general we write 

G = G a p = G a 5 a p + C a p (2.11) 

Of course, in the paramagnetic region there is no explicit replica dependence (in a 
RSB phase C a p however depends^ 13-15 ! upon the ct,/3 overlap). 

The observables calculated by direct averaging as in 2.1, are recovered via 

G = G a \ n — >0 

C=C a(3 \ n _^ Q (2.12) 
Under RSB, one can relate^ 14 ' 15 ! them by, 



r=Ti E c ° 



C — , -T- / u a/ 3 
n(n — 1 



(2.13) 

n >0 



3. LEGENDRE TRANSFORM 

To keep things simple we work from the paramagnetic phase i.e. with (<f a ) n = 0. 
We can then forget about the first Legendre transform that takes H into (<fi a ) n and 
concentrate on the second transform! 21-23 ] taking A a p(p) into ((p a (p)(pp(—p)) . We 
treat here A as a source which we extend to values A = A a p(p) for the purpose of 
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generating appropriate observables, with in the end A aj g(p) — > A. Together with 
W n = In ( Z n ^ we introduce 

W n {A} = -r n {G} + i tr A G (3.1) 



with 



and 



GaM = , A 7 , = G a (p)<W + C a/9 (p) (3.2) 



instead of (1.4,5). 

The T n functional is itself given by 



-r n {G} = 1 tr In G - 1 tr 1 G + /C«{G} (3.4) 



that is exhibiting components, 



1 1 

-T n {G a ;C a p} = -^lnG a + - tr In (<5 Q/3 + G" 1 ^) 

- lYl i G °] _1 + c ««\ + E ^ c^} ( 3 - 5 ) 



a s = l 

Here /C 1 -* is the 1-irreducible functional built with vertex and G^p lines (i.e. such 
that by cutting off two such lines, whether connected (G a 5 a p) or field connected (C a p) 
the representative graph does not fall into two disconnected pieces). The subscript s 

in fCs is the number of free replica indices, after account of the 5 a p constraints of 
the connected propagators. 



3.1 Stationarity condition 

We consider separately, stationarity with respect to off diagonal and diagonal 
components. 

Off diagonal component 5/5C a/ 3 : 

[l + G-'C^Gj' + A^+Eg^O (3.6) 

In the paramagnetic phase, C^p — ► C and, in the above equation, only s = 2 
contributes 

G^CG- 1 = A + -j^-{G; C} (3.7) 
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Diagonal component S/SG a = S/C aa : 

The equation obtained is more subtle to interpret because it contains both 
connected and field- connected graphs and hence provides two equations. In Appendix 
A it is shown that one equation is the Dyson equation for G a 



s=l 



skP 



= 



(3.8) 



where, in the paramagnetic phase, only s = 1 contributes. The other equation is the 
corresponding equation for C aa 



[t + G-'C] 'G^C] G^ + Aaa + Y, 



s=2 



8K 



(i) 



6G f 







(3.9) 



/—conn 



Both eqs.(3.6) and (3.9) can then be rewritten as (Appendix A) 



[l + G^C] 1 G~ 1 CG- 1 + A Q/3 + V 



SJC 



(1) 



s=2 



5C, 



a/3 



(3.10) 



an equation valid for a ^ (3 and a = (3. This seemingly formal result has the 
consequence that, contrary to what happens in the standard SG 



C(x=l- e) 







C(l) 



(3.11) 



showing that there is no jump in a RSB phase as the a D j3 overlap x is taken to be 
exactly equal to one*. 



3.2 Second-derivative matrix 

We have 

M ^^ = ecS&f) (3 ' 12) 

a matrix in p, p' and in replica pairs a[3, ^6. The structure in replica pair space has 
been analyzed by de Almeida and Thoulesst 24 5 (for the paramagnetic region and in 
the absence of diagonal components of G). Here again the dangerous sector is the 
replicon one with the matrix 

X R (p-p') = Mi - 2M 2 +M 3 = M R (p;p') (3.13) 

where 

Mi = M a/3;a/3 {p\p') 

M 2 = -Ma/3;a 7 (p]p') = M a p ]lf 3 {p\p') 

M 3 = M af3 ;- ( s(p;p') (3.14) 



* Note that for consistency, the extension A — > A a ^, introduced here, has to 
satisfy a relationship analog to (3.11). 
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The above expression simplifies greatly if one recognizes compensations occurring 
between the M components. These compensations are handily taken care of as follows. 
Consider the first functional derivative 

A., = §^ (3.15) 

A5T n ST n . . 

aa = = -T7T- ( 3 - 16 ) 

Contributing graphs are such that the ends a, (3 in A a p are necessarily field- connected, 
whereas in A aa the ends may be connected (contributing to the equation for G~ x ) 
or field-connected (contributing to the equation for C aa ). Upon a second derivative, 
consider the connectedness of the new end points to the pair of initial end points. 
These may be connected or field- connected. The structure of eigenvalues (3.11-14) is 
such that, one recovers 

X R (p;p') = M!| conn = M a p. a p(p;p% OQQ (3.17) 

where the index conn, stands for the connectedness between the right and left pairs. 

All the graphs with field- connexions between the left and right pairs compensate 
each other to only leave (3.17). From the explicit form of T n one gets, in exact form 



5 2 JCi 1] 



X R (p;p>) = G-\p)G-^ (p>) <W;0 - g SCaf)(p)5 c af){j/) 



(3.18) 



where only s = 2 contributes in the n — > limit. 

Note that X R starts with an attractive coupling making it a candidate to come 
out with a null eigenvalue. 



4. A-SUSCEPTIBILITY AND THE SG TRANSITION 

We are now in a situation that bears some analogy with the SG in field. In the 
paramagnetic region we have C^p = C, the analog (now space dependent) of the SG 
order parameter q a fs = q. As one crosses the line, defined by the vanishing of the 
lowest eigenvalue of A# (p;p') , playing the role of the Almeida-Thouless line, to avoid 
negative eigenvalues one has to break replica-symmetry! 13-15 ^ and write C a p = C(x) 
where x = a fl f3 is the overlap of the replica pair (in the Parish 16 ' sense). 

Just like in the pure system the vanishing of the jacobian signals the occurrence 
of a singularity in the insusceptibility, here the vanishing of the lowest Xr eigenvalue 
signals a singularity in the A-susceptibility. 

Ca)W ( pi p) = 5s _ sgs _^ 5 (4.1) 

Indeed in its replicon sector, we have 



Qr (p;p') = Qi (p;p') - 2£ 2 (p;p') + G s (p;p') 



(4.2) 
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which is just the inverse of Mr (p;p f ) '■ 

Gr (p;p') = [M r ]- 1 (p;p') 

that is 



Gr (p;p') = g 2 ( p ) 



Sp+p'-O + ^Mal3;al3 ip\p") 



G r (p";p') 



One may also directly write out the standard SG susceptibility 



(4.3) 



(4.4) 



XSG (ri - r 2 ) = [{<p (ri) </? (r 2 )) - (y? (ri)) (v? (r 2 ))] 2 
which is related to Gr by 



(4.5) 



(4.6) 



1,2 



It is striking to see that in the random field system, SG singularities are confined 
to 2-replica contributions i.e. looking back at (2.6) into — F 2 = W 2 — {^^j the 

free energy fluctuation and not the free energy itself. 

However, as soon as we are in a RSB phase, free energy fluctuations F2,F^... 
are no longer of order n 2 ,n s respectively but all become proportional to n and thus 
contribute in a finite way to the free energy. 

Let us see that effect on a simple example. Let us consider the lowest order 



contribution to /C 



We have 



(i) 



1,2 a,/3 



(4.7) 



a,/3 a 



>4 

J OLOL 



= n 



1-e 

dx C 4 (x) + C 4 (l) 



-/ 

Jo 



x dx — C (x) 
dx 



(4.8) 



where one has used (3.11). Hence the O (n 2 ) term is now 0(n) and one understands 
the vanishing of that contribution in the RS limit with the vanishing of the derivative 
d C(x)/dx. In general the total number of derivatives (with respect to overlaps x, y...) 
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is equal to the number of replicas involves minus one (s-lin /G lj ).This is an unusual 
example where the topology of the graphs contributing to the free energy strongly 
depends on the phase one is into. 



5. PHASE DIAGRAM 

To investigate, in the plane (A, T) what line is defined by the occurrence of a null 
eigenvalue, we first consider the Curie line. On that line the progagators are massless 
and we write them as follows, 

1 /A^\ 

C(P) = -^C _ . (5.1) 



p 2- v y p e ^ 

Here we have 9 = 2 — (fj — n) , 10 > {to = 1 in the mean field limit), and c(y) = 
if y = (the pure Ising case), c(y) = y if y — ► oo to recover a behavior in p~ 4+f i. 
Hence we may take 

A" A 

C(p) = = -f-= ■ (5.2) 

w p4- v p4—n v > 

As for the connected propagator we take first, for simplicity 

G(P) = ^ (5-3) 

noting that in the crossover region to the pure system, (5.3) will have to be modified. 

(i) Lowest order in A : on the Curie line 

To lowest order the eigenvalue equation reads, 

P 4 - 2 Va (p) - (9 A) 2 J -^oC 2 (q)h (P~ 0) = A/a (p) (5.4) 
with f\ (p) the eigenvector with A eigenvalue and 

C2{q) 



(2tt) £ V-*? g) 4 -*? q 8 ~ D ~ 



27? 



i „ ^ , s r 2 ((^-4 + ??)/2) . . 

Here C^q) can also be interpreted as the first cumulant contribution of a random 
temperature term. 

To overcome the difficulty of solving the above integral equation (or in Fourier 
transform, the "Schrodinger" equation with an "almost quartic" kinetic term) we 
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resort to the Rayleigh-Ritz variational approach that provides an upper bound, by 
writing 

A* = A = y ^p^\f(p)\ 2 -( 9 A) 2 J ^C 2 { q )m (5-6) 



f* (P)f(p-q) 



(5.7) 



Here / (p) a normalized trial wave function (i.e. (/>(0) = 1) whose parameters are 
to be determined variationally. 

Since we are looking for the lowest eigenvalue ("zero-energy bound state") we 
take / (p) = f(p) and real. We can now scale out the unit length R. 



X = 



a 



(gA) 2 c f d D q 



^4-277 R 2(D-4+fj) 

a / t n 2 b 



f d"q _J_. ( v 
J (2^ ^=D=2H^ 



- m 



R 4-2 V V^— ) R 2(D-4+ri) 

Writing the stationarity condition with respect to R, one gets 



(5. 



= 



(2-r))a {D-A + r ] ){gK) 2 b 



{R 2 ) 2 - 71 {R?) D ~ A+f) 
Hence solving from (5.8,9), we obtain finally 



(5.9) 



X=(D-6 + f] + r]) 



(D-4 + T]) (R 2 ) 2 ~ v 
with the length scale (i.e. correlation length) 



= (D - 6 + fj + i]) 



(2 - V ) (R 2 f- 4+ ^ 
(5.10) 



1 



D-4 + fjb 
2 — T] a 



x \ 2 



( 9 A) 



- D — rj — r) 



(5.11) 



As we rest on the Curie line we see that the eigenvalue upper bound remains negative 
for all A (except A* = where A = 0, a limit upon which we return below). 

It follows that, within the interval of dimension where the (ultraviolet) cut off 
does not spoil the length scaling, the paramagnetic to ferromagnetic transition is 
superceded by a paramagnetic to SG transition, provided it is meaningful to keep 
only the lowest order contribution to the Mr ip\p') kernel. 

The boundaries of the dimension interval are 



D u = 6 - f] (D u ) - fj (D u ) = 6 
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and 

D t = 4 - fj (D e ) 



Given that the standard dimensional reduction (and the associated rj = fj result) 
is no more applicable, one is entitled to take fj = 2rj which is correct^ 25-27 ! near D = 2. 
The lower critical dimension Di = 2 then obtained by using 9 = 2 — (fj — rj) in the 
vicinity of D = 2, that is rj = 1 for D = 2. 

Thus, modulo the (inessential) changes that will be introduced below for a 
treatment of the cross over region, what this very simple calculation is telling 
us reduces to the following: the results obtained by perturbation to all orders 
(dimensional reduction with 9 = 2, and fj = rj) are superceded by the occurrence 
of the SG transition which originates in the attraction existing between pairs of 
distinct replicas. In contradistinction, and a contrario, for "animals" (i.e. branched 
polymers) whose effective Lagrangean is alike the RFIM one but with a pure imaginary 
coupling! 28-30 ! ; the attraction becomes a repulsion, and in a random field, dimensional 
reduction is indeed correct t 29 l 



(ii) To test the robustness of the above result, one may follow Mezard and 
Young! 14 ! in adopting Bray's^ 17 ! approach i.e. use a screened interaction for an m- 
component system and work consistenly to a given ^ order. 

Eq.(5.6) is now replaced by 



{ghf f d D , 



m 



d D p 



-i 2 



(2n) 



nfiP) 



(p + q) 4 "^ 



Ol-o) (5-11) 



with C2(q), 4>(q) as of (5.5,7) and 



S(q) = jf-D-n-n {1 _ @ {mq _ q)] + V V e {mq _ q) 



G — D—ri—ri 
% 



(4n) D / 



T ((6-D-r,-fj) 12) 



(5.12) 

r((£>- 1 f 7y)/2)r((D-2 + Ty)/2) 

r(u-(S=p))r(V)r(^) 

(5.13) 



Here ji can either be unity, or chosen to take the best account of screening. Again 
scaling out the unit length yields 



T \2 



A 



a 



( 9 A) 



1 



R 4 ~ 2 ^ m R^(D-4+v) | J Q ( 2 7t) d \q R 



d u q 



q 



2(6-13-77-77) 



A(q) + 



r°° aD 

/ a g 2(6-13-77-77)13 



A(q) - O 



(5.14) 
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A(q) = C 2 (q)ct>(q) + 



f d D p 1 
J (2n)° JW ( p + q )4-^ 



(5.15) 



Note that the i?-derivative with respect to integration boundaries does not contribute. 
The stationarity condition upon R yields then 



o = (2- v ) 



(9A) 2 



1 



R 4-2 V 



m R2(D-4+fj) 



MoR d D g / q 



(2tt 



q R 



2(6-D-fj-r)) 



(gA) 2 (D-4 + f]) 



m 



R2 (D-4+fi) 



jiq R 



d D q 



fj, 



2(6-D-fj- 



v) A(q) . 



(5.16) 



leading again to the Rayleigh-Ritz approximation of the eigenvalue, 

! r 00 d D q 



A = 



(D- 6 + 77 + 77) {gA) 2 



2-77 



m R2(D-4+v) 



lxq R 



(2tt 



This expression remains negative for D < 6, confirming the result obtained in (i). 

(in) So far we have shown instability along the Curie line when keeping the terms 
in ^ ((7 scr A) 2 . One would obtain analogous qualitative behavior for ^ ((7 scr A) 



terms. The first repulsive contribution only occurs as S^A 



(fir scr A)^ 



-1 2 



(iv) Lowest order in A : cross over region 

If we want to use G(p) throughout the crossover to the pure limit A = 0, one 
should replace (5.3) by 



G(p) = 



p 2 v " \ p J 

with g(x) ~ x' 7_?7p as x — > and g(x) ~(7asx — ► 00, the subscript p referring to 
the pure limit, or alternatively, by 



-9 



(5.18) 



G( P ) 



P 



,2-r/p ■ 



Ai 



p 



(5.19) 



with g(0) = 1 and g(x) ~ l/x 71 Vp as x 
Now Eq.(5.8) is replaced by 



00. 



A = 



Af(i-ip) f d D p 



R 4-2 V 

and with the stationarity condition 



/ w? p4 ~ 2vf2(p) g ~ 2 ( RAi/p ) 



A 2 b 

R2 (D-4,+n) 



(5.20) 



o 



2-77 A%^- 



o) 



D — A + fj R A ~ 2r i 



dp 
(2^p 



p a ~ 2 V 2 {p) g- 2 (y/ P ) 



1 + 



2 y g(y/p) 



2 -vpg(y/p) 



A 2 b 



R 2(D-4+fj) 
(5.21) 



I 
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where y = RA^ , thus yielding, with obvious notations, 
D - 6 + fj + 7] Avb-iJ 



A = 



D - 4 + fj R^~ 2 ^ 



(g- 2 ) + ( -gg- 3 

6 — D — 7] — r) \p 



(5.22) 



A sufficient condition to keep A negative is to have g(x) be a monotonously increasing 
function, a very natural property given the above limiting values for x = 0, x — > oo. 

The correlation length R is now given by (5.21) and one verifies that when D is 
between and D u , A vanishes with A. 

(v) Lowest order in A : away from the Curie line 

As one is departing from the Curie line, the propagators become massive. In 
Eq.(5.6) e.g., the "kinetic" contribution will be increased slightly but the "potential" 
one will be sharply decreased, the mass playing the role of an infrared cutoff. Hence 
the eigenvalue will increase and become null at some point, on the SG transition line. 
To obtain that line we need the scaling functions for the propagators which are now 
significantly more complex, since they depend upon two variables x = A~® /p = Ao/p 
and y = 5T Up jp where 5T = |T — T c (A) | is the distance, for a given A to the Curie 
line. If one is willing to become more speculative, one may use a typical scaling form, 
restituting the appropriate behavior in all limits, as 

-A— v 



c -i (v) - p_i i (gr 

' } (5 ' 23) 



where 7 = v (4 — 77) , and 

G- 1 (p)=p 2 -^+p 2 -^A v - r < p + 



5T^ 



5T 7-Tp + A (T-7 P )/^ 



= p 2 -VpS 1 + x v-v p+ y 2-v P J^MJ_^!^_\ (5.24) 

Leaving out a complicated discussion to be dealt with separately, let us just consider 
the vicinity of the upper critical dimension. 

Note first that the Curie line becomes now the locus of a SG/Ferro-SG transition. 
However the G propagator being weakly dependent upon the SG order parameter, we 
shall assume it unchanged near D = 6, that is given by 

T — T c + a A ~ (5.25) 

where a is positive. 

Let us compute the SG transition line for D = 6 — e, at vanishing values of 5T 
and A. Proceeding in the same manner as in (i) and (iv) we get to leading order 

ST ~ 6A 4 / £ (5.26) 
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where b is positive and vanishes with e. We thus have 

T - T c ~ -aA + 6A 4/e (5.27) 

i.e. the SG transition line starts tangent to and remains very close to the Curie line 
for small e. 

As D decreases, the SG domain gets wider but too little is known about the 
behaviour of the Curie line itself to decide whether there is reentrance (i.e. whether 
the A exponent of the b term in (5.27) can become smaller than the one of the a 
term) . 

6. CONCLUSION 

We have shown that for small enough values of gA, one obtains a negative upper 
bound for the eigenvalues of the inverse spin glass susceptibility, as it occurs in the 
free-energy fluctuation. 

This enforces the occurrence of a SG phase and the replica symmetry breaking 
investigated by de Almeida & Bruinsma, Mezard & Young, Mezard & Monasson. 
It enforces it for all D between the upper (D u = 6) and lower (Di = 2) dimensions. 
Thus the general properties obtained via perturbation to all orders^ 2-4 !, for D < 6 
(i.e. dimensional reduction with 9 = 2) have no domain of application on the Curie 
line, if the above results are not reversed for higher values of <7 S crA. Indeed, they are 
then superceded by the SG transition for any D, 2 < D < 6. 

Whether this entails that the appropriate description is via a SG order parameter 
as in Mezard- Young remains to be seen 

(i) although the above result is unlikely to be reversed, one would like to render 
foolproof the above derivation (by extending it to all orders in g A) . 

(ii) this being given, it would also be desirable to see whether the above SG 
transition is not superceded itself by 3-replicas "zero energy bound states" , 4- replicas 
etc... (as contrasted with the 2-replica studied above). 

Or put another way, just as we have seen above, that the SG transition associated 
to the order parameter (ipafp) can be interpreted as governed by an effective random 
temperature (i.e. mass, with a 6*2(5) ~ \/q^~ D ~ 2 ^ correlation), likewise one may 
ask whether higher order parameters e.g. (ipaippip^) etc... associated with higher 
random couplings ^ will not become relevant and spoil the above results. 
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APPENDIX 



Here we show that the stationarity condition on the T n functional (3.4,5) yields 
eqs. (3.8-10) with the implication (3.11). 

Consider the contributions associated with a given graph to /CW{G}. 

Then take a given choice for the G components (with a constraint as in G a 5 a/ 3 
or no constraint as in C a p). The contributions of IC^ are then associated with 

£(1) ^(i) 

Consider now the first functional derivative 



5G n 



(Al) 



and the connectedness with respect to the input-output pair of lines (opened by 
the functional derivation). Those two lines can be connected or field- connected (if 
suppressing the A-coalescence the input-output lines fall apart) and one can always 
write 



SG n 



5Gry 



SG n 



(A.2) 



field - 



With (A.2) one can separate out the stationarity condition with respect to G a (or 
C aa ) as (3.8) and (3.9). 

By inspection, one then writes 



5G a 



6JC 



(i) 

s+l 



field— conn a—f3 

Noting that any 5/5C a p contribution is by definition field-connected. Hence one 
obtains (3.10) and the result that the equation for C aa (3.9) is obtained by working 
with a pair of distinct replicas, and letting {after functional derivation) j3 — >■ a, 
implying (3.11). 
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